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ABSTRACT 

We investigate the critical core mass and the envelope growth timescale, as- 
suming grain-free envelopes, to examine how small cores are allowed to form gas 
giants in the framework of the core accretion model. This is motivated by a the- 
oretical dilemma concerning Jupiter formation: Modelings of Jupiter's interior 
suggest that it contains a small core of < IOM0, while many core accretion mod- 
els of Jupiter formation require a large core of > lOM® to finish its formation 
by the time of disk dissipation. Reduction of opacity in the accreting envelope is 
known to hasten gas giant formation. Almost all the previous studies assumed 
grain-dominated opacity in the envelope. Instead, we examine cases of grain-free 
envelopes in this study. Our numerical simulations show that an isolated core of 
as small as 1.7M® is able to capture disk gas to form a gas giant on a timescale 
of million years, if the accreting envelope is grain-free; that value decreases to 
O.75M0, if the envelope is metal-free, namely, composed purely of hydrogen and 
helium. It is also shown that alkali atoms, which are known to be one of the 
dominant opacity sources near 1500K in the atmospheres of hot Jupiters, have 
little contribution to determine the critical core mass. Our results confirm that 
sedimentation and coagulation of grains in the accreting envelope is a key to 
resolve the dilemma about Jupiter formation. 

Subject headings: accretion, accretion disks — planets and satellites: formation 
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Introduction 



Core masses provide clues to unveiling the origins of gas giants. In the core accretion 
scenario for gas giant formation, a forming solid protoplanet (i.e., a proto-core) experiences 
a rapid gas capture from a pr otoplanetary disk to form a massive gas env elope, when its 



mass exceeds a critical mass (Mizuno 



1980 



Bodenheimer fc Pollack 



19861 ). That critical 



core mass must be reac hed within the lifetime of the disk gas of several million years 



(e.g. 



Haisch et al. 



200ll ). which places a limit on the core mass of a formed g as giant. 



Beca.use of the slow incr ease in the critical core mass with core accretion rate fIStevenson 



1982 



Ikomaetal 



20001 ) . faster formation in general results in larger core mass. Indeed, 



i n many core-accret i on models that a r e successful in form ing Jupiter within several Myr 



(jPoUack et al. 



1996 



Inaba et al. 



2003 



Alibert et al. 



2005 



etc.), the resultant core mass is 



as large as ~ lOM^ or more. 



In contrast, the mass of Jupiter's present core is inferred to be small. 



Saumon fc Guillot 



( 120041 ) made an extensive investigation of the interior structure of Jupiter, finding successful 
models that are consistent with the observed values of its gravitational moments and 
equatorial radius by using a variety of equations of state (EOSs) for hydrogen and helium. 
They demonstrated that the possible core mass of Jupiter is smaller than ~ 10 Mq. This 
is also sup ported by recent calcula tions with an ab initio EOS derived in the first-principle 



ap proach fjNette 



by 



Militzer et al. 



mann et al 



(2008 



20081 ). While a more massive core of > lOM® is reported 



who used their own EOS of hydrogen- h eliurn mixtures based on 



density functional molecular dynamics. 



Fortney fc Nettelmamj (j2009[ ) pointed out the 



difference in the mass fraction of helium used by the two groups is responsible for this 
discrepancy in the derived value of Jupiter's core mass. Although this pending problem 
about Jupiter's core may arise from the uncertainty of EOS, the core mass suggested by 
interior modeling is, on an average, smaller than that derived by formation theories. This 
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fact motivates us to know how small a core can start the rapid gas accretion to form a 
massive envelope within several Myr. 

Reduction of opacity in the protoplanet's envelope has the potential to make a 
small core possible. As the opacity becomes small, the critical core mass decreases and 
the post-critical- mass gas accretion becomes fast, because low opacity in the envelope 
makes it difficult to maintain the envelope's h ydrostatic structure without gravitational 



energ y released by contraction of the envelope (IMizuno 



1980 



Stevenson 



1982 



Ikoma et al 



2000|). Since the opacity sources are dust grains and gaseous components in the envelope, 
a minimum critical core mass is achieved in the case of grain- free envelopes. All the 



Mizund (119801 ) (see Section [2]) assumed 



previous studies except one calculation done by 
grain-dominated opacity in the outer envelope. Thus, in this paper, we consider grain-free 
envelopes and make an extensive investigation of the critical core mass and timescale for 
gas accretion. 



Results from the recent works by Podolak and his colleagues (jPodolakI 12003 



Movshovitz fc Podolak 



20081 ) are encouraging. They have directly simulated the dynamical 



behavior of dust grains to determine their size distribution, and then calculated grain 
opacity in the accreting envelope. Their numerical simulations revealed that grain opacities 
in the envelope can be much lower than those in the protoplanetary disk. The reason is that 
small grains initially suspended in the outer envelope quickly grow large in size and then 
settle down into the deep envelope where temperature is high enough that grains evaporate. 

In the following section, the details of the gas opacity used in this study are described. 
Our numerical calculations and results for the critical core mass and the timescale of gas 
accretion are shown in Sections 3 and 4, respectively. We discuss a possibility of Jupiter 
formation with a small core in Section 5. 
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2. Gas Opacity 



The opacity of tlie envelope gas is lowest when the envelope contains only hydrogen 
and helium. In this study, we first consider such a metal-free casqj. We compute chemical 
equilibrium between H2, H, H+, H^, H^, H^, and e^, and then calculate the opacity that 
includes the bound- free and free- free absorptions by H~, Rayleigh scatterings by H2, H, and 
He, Thomson scattering by e~, and the collision-induced absorptions (CIA) due to H2— H2, 
H 2— H, H2— He, and H— He. Values of quantities relevant to those calculations are given 



m 



Lenzuni et al 



(Il99l[ ) and references therein. The CIA opacities are computed 



the latest programs and tab 



Borvsow fc Frommhold 



1996 



Borysow et al. 



1990 



1997 



es availab 



Borvsow 



2000 



e on Borysow's web page 



1991 



Borysow 



1992 



20021 ). 



jy usi ng 



(IBorvsow et al 



Zhen fc Borysow 



1995 



1985 



Birnbaum et al. 



Metal-containing envelopes with no grains are also considered in this study. It is 
uncertain what fraction of metals in protoplanetary disks are incorporated in dust grains 
and what faction remains in their gaseous forms. For example, the amount of water 
adsorbed onto grain surfaces depends strongly on the thermal state of the protoplanetary 



disk. Water is adsorbed onto grain surfa ces in a relative 



in the gas phase in a relatively hot disk (IMarkwick et al. 



y CO. 



2002 



d disk, while it remains 



Woods fc Willacv 



20091 ). 



Adsorpti on rates of water are also sensitive to the abundances and sizes of dust grains in 



the disk ( lAikawa fc Nomura 



20061 ). but such dust properties remain poorly known. Alkali 



atoms may be important: They are known to be one of the dominant g as opacity sources 



near 1500K in the ou ter atmospheres of brown dwarfs and giant planets fiGuillot et al. 



Burrows et al. 



1994 



20001 ) ■ though alkali atoms have a tendency to reside onto grain surfaces 



^In the community of astrophysics, the word 'metal' represents elements heavier than 
helium. 

^http:/ /www. astro. ku.dk/ aborysow/programs/index.html 
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in protoplanetary disks (IHasegawa fc Herbstlll993l ). Therefore, we consider two cases in 
addition to the metal-free case: In one case called the alkali case hereafter, alkali atoms are 
all in the gas phase; in another case called the no-alkali case, they are absent in the gas. 



Freedman et al. 



(120081 ) 



In the alkali case, we use the gas-opacity data provided by 
which incorporate the revised solar abundances and alkali atoms. The solar abundances 
have been recently revised by re-analyses of two forbidden lines of neutral oxygen and 
carbon from the solar photosphere under the assumption of the local thermodynamic 
equilibrium, although their values should still require s crutiny because of uncerta i nties in 



3-D hydrodyn amical models for the solar atmosphere (lAUende Prieto et al 



Lodder 



2001 



2002 



20031 ). Opacity data without alkali atoms are kindly provided by D r. Freedman (in 



person al communication). In practice, we use the up-to-date opacity data of 



Freedman et al. 



( 120081 ). which cover the wider ranges of temperature and pressure compared to the published 
ones and include the latest HITRAN spectroscopic data and an additional opacity source, 
CO2 (in personal communication). 



Critical Core Mass 



The spherically-symmetric hydrosta tic structure o f a protoplanet is simulated in a 



manner similar to previous studies (e.g. 



Mizuno 



19801 ). The protoplanet consists of a 
solid core and a gaseous envelope. The core has a constant density of 3.2g cm~'^ and its 
structure is not computed. The envelope is assumed to be in purely hydrostatic equilibrium 
and have a uniform chemical composition. While the composition de pends on the c ase, 



the abundance ratios of elements taken into account are always solar (ILodders 



2OO3I ). All 



simulations are cor nputed with sp l ine-in terpolated SCVH EOS tables for hydrogen and 



helium provided by 



Saumon et al 



(Il995[ ) , contributions from heavy elements being ignored. 



The inner boundary conditions are applied at the core surface. At the outer boundary, the 
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temperature and density are equal to the midplane values of a protoplanetary disk at the 
protoplanet's orbit, T^^sk and Pdisk, because the envelope is assumed to be in equilibrium 
with the disk gas at the outer edge. The outer radius of the protoplanet is defined by the 
smaller of the accretion radius and the tidal radius. Note that our results are insensitive 
to the outer boundary conditions. We handle the accretion ra te of planetesim als, Mc, as a 



free parameter in this study. We follow the same proc edure as 



the critical core mass, Merit (please see 



Mizund (119801 ) to determine 



Mizund (119801 ) in detail). Table 1 summarizes input 



parameters and their values used in this study. 

Figure 1 shows M^it as a function of Mc for the three cases. For comparison, we 
plot results for additional two cases in which the envelope contains both grains and gas 
as opacity sourc es. We assume that the grain opacity takes / times values of the grain 



opacity given by 



Pollack et al. 



(119851 ): their calculations assumed a nearly interstellar size 



distribution of dust grains. 

Reduced opacity results in small M^rit, as shown in Figure 1. For instance, we now 
focus on results for M^ = 1 x lO~^M0/yr. In the cas e of f = 0.01 (double dot-dashed line). 



Hubickyj et al. 



20051), Merit = lOMe. 



which is often adopted in core accretion models (e.g. 
Removal of grains lowersM^rit to a few Mq; 3.5M^ for the alkali case (dashed line) and 
1.5M^ for the metal-free case (solid line). 

Compared to the alkali case, the metal- free case always produces smaller Merit- This is 
because molecules composed of oxygen and/or carbon such as H2O and CO2 are effective 
opacity sources as for determination of Merit- In contrast, comparison between the results 
for the alkali (dashed line) and no-alkali cases (dotted line) demonstrates that alkali atoms 
have little contribution to determine M^it- The reason is that convection governs heat 
transfer in deep, hot parts of the envelope where alkali atoms have a great contribution to 
the opacity. 
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Figure 1 also demonstrates that all the lines except one for the metal-free case are 
converging, as decreases. In the case of low Mc, namely, low luminosity, the outermost 
isothermal layer extends deep in the envelope. In that isothermal layer, the density increases 
rapidly to keep the pressure gradient needed for supporting the core's gravity. Thus, 
because of high densities, the opacity from H2O and CO2 dominates the grain opacity in 
the deep radiative envelope in spite of low temperature. 



Mizunol fll980h 



Finally it is worth mentioning differences from previous studies, 
calculated M^it in the metal-free case for only one value of core accretion rate, 
1 X lO'^M^/yr, and derived l.SM^. While the CIA opacities at high temperatures have 
been revised thanks to progresses of quantum mechanical models and improvement of 
experimental data in the 1990s, Miz uno's (1980) value of M^nt is in good agreement with 



our result for the same parameters. 



IkomaetaL 



(I2OOOI ) investigated cases with low grain 



opacity. Ou r calcu lations for / = 0.01 and / = 0.001 yield larger Ma-it than those from 



IkomaetaL 



(l2000l ). because they assumed relatively small grain opacities compared to 



those used in this study. 



4. Envelope Growth Timescale 



The critical core mass decreases as the core accretion rate decreases, as shown in Fig.l. 



There 



:ore, the case of Mr, — )• (i.e., isotherm al envelope) yields the absolute minimum of 



Merit (I Sasaki 



1989 



Pecnik fc Wuchterl 



for the core to capture disk gas (jlkoma et al. 



20051). However, for smaller Mrrit. it takes longer 



2000 



Ikoma fc Genda 



20061). Because there 



20011 ) 



are time constraints on gas giant formation such as the disk's lifetime ( iHaisch et al 
a practical minimum of M^nt is determined in this respect. It is necessary to estimate how 
long it takes for a protoplanet with a given core mass, Mcore, to capture disk gas. 
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We consider an isolated protoplanet, namely, accumulation of the envelope after 
planetesimal accretion is halted. The isolated protoplanet always experiences contraction 
of its envelope and captures disk gas, because of no energy supply due to planetesimal 
accretion. We simulate the quasi-static evolution of the envelope with a given M^ore after 
planetesimal accretion is halted and evaluate the growth timescale of the envelope. As the 
growth timescale of the envelo p e, we present values of the characteristic growth time, Tg, 



that is defined by 



Ikoma et al. 



mm . 



Figure 2 plots Tg as functions of Mcovc- For comparison, the results for / = 0.01 
and 0.001 are also shown. The growth time Tg is found to be much shorter in the case 
of the grain-free opacities than Tg for the grain-dominated opacities. For example, when 
M(. = 3M0, Tg = 7 X 10^ yr in the metal-free case (solid line) and 7 x 10^ yr in the alkali case 
(dashed line), while 2 x 10^ yr for / = 0.01 (double dot-dashed line). The envelope growth 
is regul ated by the Kelvin -Helmholtz contraction of the envelope; that is, Tg is proportional 



to 1/L (jlkoma et al. 



20001), where L is the luminosity at the envelope's outer edge. Since a 
change in opacity compensates that in luminosity, k being an averaged opacity, Tg should 
be proportional to R/L. In the cas e of the grain-domi nate d opacities, the d e pende nce is 



simple, namely, Tg cx /, as shown in 



Ikoma et al. 



mm and 



Ikoma k. Gendal (120061 ). Figure 



2, however, demonstrates that such a simple scaling is inappropriate in low R cases. The 
curves for grain-free cases are steeper than those for grain-dominated cases; the reason has 
been already described in the previous section. 

In any case, we have found that the envelope growth time is significant shorter in the 
case of the grain-free opacities. Based on this fact, we discuss the minimum critical core 
mass from the viewpoint of gas giant formation, and the possibility of Jupiter formation 
with a small core in the following section. 
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5. Discussion 



There is a theoretical dilemma concerning Jupiter formation, as described in 
Introduction. Modelings of Jupiter's in terior suggest that Jupiter has a small core of 



< lOMe (e.g., 



Saumon fc Guillot 



2004 ). while many core-accretion models of Jupiter 



formation r equire a large core of > lOMm to finish its 



dissipation (jPoUack et al. 



1996 



Alibert et al. 



2004 



2005 



brmation by t 



le time of disk 



Fortier et al. 



2007 



20091). The 



disk instability scenario has 



Boss 



2000 



Mayer et al. 



)een revisited as an alternative scenario of their formation (e.g. 



2002) 



In this study, we have demonstrated that reduced opacities in the protoplanet's envelope 
have the potential to resolve this dilemma. From Fig. [21 one finds that M^orc = O.TSM^ for 



Tg = 1 Myr in the metal-free case; 



1.7Mffi even in the alkali case. Given observed 



lifetimes of protoplanetary disks of several Myr, the fact above indicates the reduction of 
opacity allows Jupiter to have a small core that is consistent with interior modelings in 
principle. The feasibility of such minimum M^it depends on opacities in the protoplanet's 
envelope, while this does not change our conclusion that the minimum Merit obtained here 
provides the lowest limit to core masses of gas giants to which the core accretion model can 
apply. We need more extensive investigation of gas giant formation into which sedimentation 



and coagulation of grains in the accreting enve 



opacity was already pronounced (jPodolak 



2003 



ope are incorporated, a. 



Movshovitz fc Podolak 



though reduction of 



20081 ). 



Our results also shed light upon growth of s olid cores. Many co re accretion models 



Pollack et al. 



19961 ). However, a gas giant 



assumed the presence of a single protoplanet (e.g 
is in practice thought of as being formed in a system of multiple protoplanets embedded in 
a protoplanetary disk. Compared to cases of a single protoplanet, the final mass of a core 



should be sm all in the case of a multiple-protoplanet system. According to 



Kokubo fc Ida 



(11998 



2000l ). the isolation mass is a few around 5AU. Even such a small core is enough 



for a gas giant to capture disk gas within several Myr, as demonstrated in this study. We 
need to perform comprehensive simulations on gas giant formation in multiple-protoplanet 
systems, which incorporate adequate models of pla netary accretion such as fragmentation of 



planetesimals and e-damping due to gas drag (e.g. 
will be our future work. 



Duncan et al. 



20091 ). These calculations 
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Fig. 1. — Critical core masses, Mcnt, as functions of core accretion rate, M^. The solid, 
dashed, and dotted lines represent the results for the metal-free, alkali, and no-alkali cases 
(see Section [2] for the definitions). For comparison, the results of two cases in which there 
exist grains in the envelope are also shown with the double dot-dashed and dot-dashed lines: 



/ = 0.01 and 0.001, respectively, where / is the grain depletion f actor (see the text 



definition). We have used the grain opacity tables derived from 



Pollack et al. 



or its 



(119851). [A 



color version of this plot is available in the electronic edition of The Astrophysical Journal] 
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Mcore (Me) 

Fig. 2. — Typical growth timescale of the envelope, Tg (see the text for the definition), 
as a function of Mcore- As in Fig.l, the solid, dashed, dotted lines represent the results 
of the metal-free, alkali, and no-alkali cases respectively. The results of / = 0.01 (double 
dot-dashed line) and 0.001 (dot-dashed line) are also shown. [A color version of this plot is 
available in the electronic edition of The Astrophysical Journal] 
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Table 1: Input Parameters and Their Values. 



Parameter 

semimajor axis, a 

core density 

disk temperature, T^isk 

disk density, pdisk 

planetesimal accretion rate, Mc 



Value 
5.2 AU 
3.2 g/cm^ 
150 K 

5.0 X 10"^^ g/cm^ 

1.0 X 10^3 to i Q X 10-10 Me/yr 



